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)Cl ■ Abstract 

The model-independent radiative corrections to deep-inelastic scattering of unpolar- 
^ ■ ized electron beam off the tensor polarized deuteron target have been considered. The 

CN . contribution to the radiative corrections due to the hard-photon emission from the elastic 

I electron-deuteron scattering (the so-called elastic radiative tail) is also investigated. The 

^T) ' calculation is based on the covariant parametrization of the deuteron quadrupole polar- 

CD ■ ization tensor. The numerical estimates of the radiative corrections to the polarization 

observables have been done for the kinematical conditions of the current experiment at 

. HERA. 
^ ■ 

j^' PACS: 12.20.-m, 13.40.-f, 13.60.-Hb, 13.88.+e 



1 Introduction 



The flavor structure of nucleons is described in terms of parton distribution functions. Most 
of the information on these functions has up to now come from inclusive deep-inelastic scatter- 
ing process: experiments where only the scattered lepton is detected. The investigation of the 
nucleon spin structure involves now new types of reactions. For example, the HERMES exper- 
iment was specifically designed to perform accurate measurements of semi-inclusive reactions, 
where besides scattered lepton also some of the produced hadrons are detected p. 

The polarized nuclei of deuterium and helium-3 are used to extract information on the 
neutron spin-dependent structure function gi{x) 0. In analyzing the experimental data on 
inclusive spin asymmetries for deuterium in order to deduce the spin-dependent structure 
function gf one should to take into account a small effect due to possible tensor polarization in 
this spin-one target. This is connected with the presence in a deuteron target of an additional 
tensor polarized structure functions [1^ . So far, the spin-structure studies have been focused on 
the spin-1/2 nucleon. Different spin physics exists for higher-spin hadrons such as the tensor 
structure in the deuteron. The measurement of these additional spin-dependent structure 
functions provides important information about non-nucleonic components in spin-one nuclei 
and tensor structures on the quark-parton level jn|. A general formalism of the deep inelastic 
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electron-deuteron scattering was discussed in Ref. where new four tensor structure functions 
bi{x), i = 1 — 4 were introduced. They can be measured using tensor polarized target and 
unpolarized electron beam. Among these new structure functions, only one, bi, is leading twist 
in QCD m and it was found that this function is small for a weakly bound system of nucleons 
(for example, deuteron). Therefore, the measurement of hi for the case of the deuteron can 
give information about its possible exotic components. 

From the theoretical point of view the spin-dependent structure function bi{x) was inves- 
tigated in a number of papers. The available fixed targets with J > 1 are only nuclei (the 
deuteron is the most commonly used nucleus). If the nucleons in the deuteron are in an S state 
then bi{x) = 0. For nucleons in a D state, bi{x) 7^ in general It was found [Sj that in 
a quark-parton model the sum rule / dxbi{x) = is generally true if the sea of quarks and 
antiquarks is unpolarized (and it was shown how this sum rule is modified in the presence of a 
polarized sea). Mankiewicz IHj has studied the bi{x) for the p meson and noticed empirically 
that / dxbi{x) = in his model. It was shown in Ref. that multiple scattering terms at 
low X can still lead to 61 7^ even if only the S'-wave component is present. Various twist 
two structure functions of the deuteron (in particular, bi) have been calculated in a version 
of the convolution model that incorporates relativistic and binding energy corrections [8^. A 
simple parametrizations of these structure functions are given in terms of a few deuteron wave 
function parameters and the free nucleon structure functions. The tensor structure functions 
were discussed in Ref. P for the case of lepton scattering and in hadron reactions such as the 
polarized proton-deuteron Drell-Yan process. 

As it is known, the HERMES experiment has been designed to measure the nucleon spin- 
dependent structure functions from deep-inelastic scattering of longitudinally polarized positrons 
and electrons from polarized gaseous targets {H, D,^ He). In 2000, HERMES collected a data 
set with a tensor polarized deuterium target for the purpose of making a first measurement of 
tensor structure function bi{x). The preliminary results on this structure function are presented 
in Ref. HU] for the kinematic range 0.002 <x < 0.85 and O.lGeV^ < < 20GeV^. The pre- 
liminary result for the tensor asymmetry is sufficiently small to produce an effect of more than 
1% on the measurement of gf. The dependence of bi on x variable is in qualitative agreement 
with expectations based on coherent double scattering models [TT[ [Tallin] and favors a sizeable 
value of bi at low x region. This suggests a significant tensor polarization of the sea-quarks, 
violating the Close-Kumano sum rule j5|. 

The radiative corrections to deep-inelastic scattering of unpolarized and longitudinally po- 
larized electron beam on polarized deuteron target were considered in Ref. for the par- 
ticular case of the deuteron polarization (which can be obtained from the general covariant 
spin-density matrix ^3] when spin functions are the eigenvectors of the spin projection op- 
erator). The leading-log model-independent radiative corrections in deep-inelastic scattering 
of unpolarized electron beam off the tensor polarized deuteron target have been considered in 
Ref. [16j. The calculation is based on the covariant parametrization of the deuteron quadrupole 
polarization tensor and use of the Drell-Yan like representation in electrodynamics. 

Current experiments at modern accelerators reached a new level of precision and this cir- 
cumstance requires a new approach to data analysis and inclusion of all possible systematic 
uncertainties. One of the important source of such uncertainties is the electromagnetic radia- 
tive effects caused by physical processes which take place in higher orders of the perturbation 
theory with respect to the electromagnetic interaction. In present paper we give the covariant 
description of deep-inelastic scattering of unpolarized electron beam off the tensor polarized 
deuteron target (the polarization state of the target is described by the spin-density matrix of 
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the general form) taking into account the radiative corrections 

e~{h)+d{p)^e-{h)+X{p,). (1) 

The corresponding approach is based on the covariant parametrization of the deuteron 
quadrupole polarization tensor in terms of the 4-momenta of the particles in process (1) [16j . 
We performed also the numerical calculations of the radiative corrections for the kinemati- 
cal conditions of the experiment PUJ- The contribution of the radiative tail from the elastic 
erf— scattering is considered separately. 



2 Born approximation 



The standard set of variables which is usually used for the description of deep-inelastic 
scattering process is 

X = y = V = 2pki, = -Vxy, q = h - /ca, (2) 

where q is the 4-momentum of the intermediate heavy photon that probes the deuteron struc- 
ture. To begin with, we define the deep-inelastic scattering cross section of the process (1) in 
terms of the leptonic L^^^, and hadronic IV^i, tensors contraction (in the Born approximation we 
can neglect the electron mass) 

_^^_na^y 

dxdQl ~ VQ^x^" ^'^ ■ ^ ^ 

Note that only in the Born approximation (without accounting radiative corrections) q = 
h - k2, Q\ = -q^ = 2kik2. 



The Born leptonic tensor (for the unpolarized case) is 



= <f9tiy + ^hfihu + ki„k2fj.). (4) 



The hadronic tensor is defined by the following way 

W^, = (27r)3 ^ 6'^'\k, +p-k2~ p.)j;j^, 

X 

where is the electromagnetic current for the j* + d —>■ X transition (7* is the virtual 
photon). The sum means summation over the final states and the bar means averaging over 
the polarizations of the target and summation over the polarizations of the final particles. To 
write down the hadron tensor in terms of the structure functions we define first the deuteron 
spin-density matrix (further we do not consider the effect caused by the vector polarization of 
the deuteron) 

2 P^Pu ^ 

here and Q^j, are the target-deuteron polarization 4-vector and the deuteron quadrupole- 
polarization tensor. The corresponding hadron tensor has both the polarization-independent 
and polarization- dependent parts and in general case can be written as 

W,, = W,,{0) + W,,{V) + W^,{T), (6) 
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where iy^,y(0) corresponds to the unpolarized case and Wfj_^(y){Wfj,^{T)) corresponds to the case 
of the vector (tensor-)-polarization of the deuteron target. The Wf^u{0) term has the following 
form 

here M is the deuteron mass, and Wi^2 are the unpolarized structure functions depending on 
two independent variables x and q^. The part of the hadron tensor which is dependent on the 
quadrupole polarization tensor can be represented as 

W^lu{T) = J—^iQapqaQpiBig^^ + —Pf^Pu) + B^q^ip^Qua + PuQfla) + pqBiQ^u}- (8) 

Here Bi{i = 1,2,3,4) are the spin-dependent structure functions (caused by the tensor polar- 
ization of the target). They are also functions of two variables and x. Since the hadron 
tensor W^jy(T) is symmetric under /i i/ the electron beam does not have to be polarized for 
measuring these new structure functions. 

We used the following notation in formula (8) 

P\ ^ I QfJ-Qu ^ QuQa ^ QfiQa ^ ^ PI /qn 

Note that the deuteron spin-dependent structure functions Bi are also related to the structure 
functions bi, introduced in Ref. jlj, in the following way 

B, = -b„ 52 = | + 63 + &4, 53 = |-^, S4=|-&3. (10) 
6 b 2 6 

When calculating radiative corrections it is convenient to parametrize the polarization state 
of the deuteron target in terms of the 4-momenta of the particles participating in the reaction 
under consideration. Therefore, first, we have to find the set of the axes and write them in 
covariant form in terms of the 4-momenta. If we choose, in the laboratory system of reaction 
(1), the longitudinal direction 1 along the electron beam and the transverse one t in the plane 
(ki,k2) and perpendicular to 1, then 

^(i) _ 2rfci^ - _ fc2M - (1 - y - 2x|/r)fci^ - xyp^, _ '2e^xpaPxkipk2a , . 

- M ' ~ d ' ~ Vd ' ^^^^ 



d = ^Vxyb, b = l- y- xyr, r = MyV. 

We chose one of the axes along the direction I because in the experiment on measuring the 
bi structure function ^U] the direction of the magnetic field, used for the polarization of the 
deuteron target, is along the positron beam line. The direction of the magnetic field provides 
the quantization axis for the nuclear spin in the target. 

One can verify that the set of the 4-vectors S"^''*'"-' satisfies the following properties 

Si^»)sW = -S^p^ 4% = 0' c^,P = l,t,n. (12) 
One can make sure also that in the rest frame of the deuteron (the laboratory system) 

5i'^ = (0,l), 5W = (0,t), 5W = (0,n), 
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n2-(nin2)ni ni x na ki,2 

Vl-(nin2)' Vl-(nin2) l^i'^l 

If to add one more 4-vector 5"^°^ = Pn/M to the set (11), we receive the complete set of the 
orthogonal 4-vectors with the following properties 

gMgim)^g^^^ ^M^(n)=^_^ m, TZ = 0, /, t, n. (14) 
This allows to express the deuteron quadrupole polarization tensor, in general follows 

Q^iu — Sjj^'' si"''' Rmn = Sj°'''Sl^''Rai3, Ra/B — R/Sa, Raa — (15) 

because the components i?oo, Roa and Rao identically equal to zero due to condition Q^uPv — 0. 

In the Born approximation the components Rin and Rtn do not contribute to the cross 
section (since the 4-momenta and A'l^j are orthogonal to the 4-vector S^^^) and the expansion 
(15) can be rewritten in the following standard form 

Q,^ = [Sll^Si'^ - i^W^W]/?^ + \si^^Si'\Ru - Run) + (^^^4*^ + Si^S^^)Ru , (16) 

here we took into account that Ru + Ru + Rnn = 0. 

In further, the deep-inelastic scattering of the unpolarized electron beam from the tensor 
polarized deuteron target is considered. Thus, we have to calculate only the convolution of the 
Born leptonic tensor L^^^ and hadron tensor Wni,{T) caused by the tensor polarization of the 
target 

S''{T) = L^,W^,{T) = 8^{-l[xy^B, + (a - 1 + y)B2 + yBs]Qo + (17) 

y y 

+J[(2 - y)B^ - yB4Qi + B^Qu}, 

where a = xyr, Qo = QafsQaqp, Qi = Qa/sqah/s, Qu = Qafshahp- Using the formulae for the 
vectors Sj^^ we can calculate the convolutions. After simple calculations we have 

^ ^ [SllRll + Stt(Rtt — Rnn) + SltRlt], (18) 



with 



Sii = [2xbT - y{l + 2xr)2]G + 26(1 + 3xt)B3 + {b - a)B^ 



XUT 

Sit = 2y — [2(y + 2a)G + {2 - y - Ah)B^ + yB^], (19) 
Stt = -2xbr{G + S3), G = xyBr - -B2 . 

y 

So, in the general case, the cross section of deep-inelastic scattering of unpolarized elec- 
tron beam from the tensor polarized target is determined, in the Born approximation, by the 
components of the quadrupole-polarization tensor Ru, Ru and the combination {Ru — Rnn)- 

Consider just one more, commonly used, choice of the coordinate axes: components of the 
deuteron polarization tensor are defined in the coordinate system with the axes along directions 
L, T and N in the rest frame of the deuteron, where 

^-W^V T-^^^^i^, N = n. (20) 
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The corresponding covariant form of set (20) reads 

ciD _ 2r(A;i - fca)^ - yp^ o(r) _ + 2xr)A;2M - (1 - 1/ - 2xr)A:i^ - x{2 - y)p^ 

~ ' - VVxbh ^ ' 

and the expansion of the deuteron polarization tensor is defined in full analogy with (16) 

Q.. = [Sl'^Sl^^ - ls^pSP]R,, + ^-S^;^^sP{Rtt - Rnn) + (4'^^f ^ + Spsi^^)R,T . (22) 

These two sets of the orthogonal 4- vectors are connected by means of orthogonal matrix which 
describes the rotation in the plane perpendicular to direction n = N 

^i^) = cos^^« + sin^SW, ^f) = -sin^^«+cos^4*), (23) 



y(l + 2xT) ^ xbr 

cose = — — , sin^ = -2i' 



'yh V h 

The part of the differential cross section that depends on the tensor polarization can be 
written as follows in this set of axes 



[SllRll + Stt{Rtt — Rnn) + SltRlt] , (24) 



Sll = -hG + 2hB; + ^[(1 - y){y - 2xr) - 2a{y + xr)] 



Stt = '^B,, Slt = 2^^{2 - y){Bs + ^B,) . (25) 



3 Radiative corrections 



In this work we consider only QED radiative corrections to the deep-inelastic scattering pro- 
cess (1). We confine ourselves to the calculation of the so-called model-independent radiative 
corrections when the photons are radiated from the lepton line and the vacuum polarization is 
also taken into account. The reason is that it gives the main contribution to radiative correc- 
tions due to the smallness of the electron mass, and can be calculated without any additional 
assumptions. Nevertheless, these radiative corrections depend on the shape of the deuteron 
structure functions (both spin-independent and spin-dependent) by their dependence on the x 
and variables. 

There exist two contributions for radiative corrections when we take into account the cor- 
rections of the order of a. The first one is caused by virtual and soft photon emission that 
cannot affect the kinematics of process (1). The second one arises due to the radiation of a 
hard photon 

e-ih) + d{p) ^ e-ik2)+^ik)+X{p,) . (26) 

The leptonic tensor, corresponding to the hard-photon radiation, is well known [T7\ ITHj . 
For the case of unpolarized electron beam it can be written as 

= ^oQpiv + AikiJ^i^ + A2k2t,k2u , (27) 
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^0 = — 2m g ( — + — , Ai = -4 \ ^ 



where 



X1X2 \Xi X2/ \X1X2 X2 

. J 2m^ \ r 7 / n 

VX1X2 Xl / ^ g2 A- 

with Xi 2 = 2kki 2, m is the electron mass, = X2 — Xi ~ Q%j and in this chapter q = 
ki — k2 — k. The hadronic tensor in this case has the same form as the hadronic tensor in the 
Born approximation, but momentum transfer q differs from the Born one and the structure 
functions Bi depend on the new momentum q. Here and further we neglect the terms which 
are zero at m — >■ 0. 

We consider the hard-photon (with the energy u > Ae, where A << 1) emission process 
using the approach of paper [THj where it was applied to the process of deep-inelastic scattering 
on unpolarized target. Let us introduce the variables suitable for this process 

— + 2pq . —q^ xyr ^ , , 

^ = 77 = T7 ' ^ = "TiF' ^ = = ^1.2 = 2fcfci,2, 

V V 2pq xyr + z 

where is the invariant mass of the hadron system produced in the scattering of the photon 
(with the virtuality g^) by the target. 

Note that the physical meaning of z variable is following: the quantity z shows the degree 
of deviation from the elastic process (erf — > erf). So, the value z = corresponds to the elastic 
erf-scattering threshold and the value z = ed/£i {£d is the deuteron bound energy, ei is the 
electron beam energy in the laboratory system) corresponds to the erf enp reaction threshold 
(quasi-elastic erf— scattering). 

The convolution of the leptonic and hadronic tensors may be represented as 

S^T) = AAo + BAi + CA2, (28) 

A = NQo[3B^ + —B2 + ^{B2 + 2^3 + ^4)], 
c Zxyr 



~t(B2 + + B,))] + VQ,\B, - %±#i<;(S3 + B,)] + ^< 



Vc — Y V 

-^(52 + 2i?3 + B,))] + VQ2[Bs{l -y) + ^%^<B, + ^4)] + -CQ22- 



with = At/V(?, c = z + xyr. 

The quantities Qq, Qi, Q2, Qn and Q22 are the convolutions of the deuteron quadrupole- 
polarization tensor and 4-momenta and they can be expressed in terms of the scalar products 
of the 4-momenta of the particles participating in the reaction and the set of 4- vectors 3^'^'""^ 
So, these convolutions are 

Qo = Qapqaqp = [{1 ■ qf -^{t- qf - ^(^ ■ qf]Ru + 21- qt- qRit + 2n ■ ql ■ qRin+ 
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+2n ■ qt ■ qRtn + ^[(^ " 1? " " <lf]{Rtt - Run). 
Qi = Qaf3qakif3 = {I ■ hi ■ Q - ■ kit ■ q)Ru + {I ■ kit ■ q + t ■ kil ■ q)Rit + I ■ kin ■ qRin+ 

+ t-kin- qRtn + ]f ■ kit ■ q{Ru - Rnn), (29) 

Qu = Qapki^kip = [(/ ■ kif - ^(t • kif]Ru + 21 ■ kit ■ kiRlt + ^{t ■ ki)\Ru - Rnn), 

Q2 = Qi{ki k2), Q22 = Quiki /ca), i-a = S^^^a^, i = l,t,n. 

Here we used the following conditions: Ru + Ru + Rnn = 0, n ■ ki = n ■ k2 = 0. For the set of 
the 4-vectors Sj^'*'""^ we have also t ■ ki = 0. 

It is convenient to separate the poles in the term {xiX2)^^ using the relation 

1 11,1 1 . 



X1X2 Qs'^-r Xi X2 

Then the radiative correction (caused by the hard-photon emission) to the differential cross 
section of deep-inelastic scattering of unpolarized electron beam by the tensor polarized target 
has the following form 

da'^ ay r d?k ^ , , , , 

' S(2;,r), (30) 



dxdQ% Vx J 2txlo 
where uj is the energy of the hard photon and 



S(^, r) = ^^{R,{z, r) + (- - -)Ri{z, r) + "^Rm^iz, r) + ^R2m{z, r)}, (31) 
Qb Xi X2 Xi Xi 



Ro = -^A, Ri = J—[(l+L)QlA--{B+C)], Rim = 2{^A-^C), i^2n^ = 2(^A-4 
It is convenient to write down the integral in Eq. (30) as 



d'^k 

S(,2, r) = + /2m + Ir, (32) 

Ztxuj 

where we separate the contributions proportional to im? 

[ d^k iq^) m? ^ , , ^ r d^k iq^) rri^ 

hm= / ^i?l„^(^,r), l2m= / 7^ 7^ :^R2m{z,r). (33) 

J 2nu Q% xi J 27raj Q% xi 

Let us consider first the integrals Jjm, i = 1, 2. In this case the numerator of the integrands in 
hmihm) is calculated in the approximation xi = 0(x2 = 0) [1^]. The hard-photon phase-space 
is written as 

d^k dz u'^dVLu , ^ , ^ 

7^ = z^=y{l~x). 34 

zTTo; z+ — z 2tt 

Using the invariance of uj'^dQk,^^ can do the integration over the angular variables dflk in 
the most suitable coordinate system, namely: in the coordinate frame ki — k2 + p = (c.m.s. 
of the scattered electron and produced hadronic system). We obtain 

uj'^dVtk 1 



2vr xh 2' 
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Integral We calculate the integrand in the approximation Xi = (besides the denom- 
inator). This approximation corresponds to the emission of the coUinear photon along the 
initial-electron momentum. In this case the variables acquire the following values 

+ ^ 2 ^2 / ^y^i 



I — xy z + xyri 

After integrating over the hard-photon angular variables, the integral Ii^ can be represented 
as follows 

hm = 7^ r -^^alNiEiiz), (35) 



Q% Jo z+- z~ 

^l{z) = T^iiiRii + T,iitRlt + T,itt{Rtt — Rnn)i 



\tt = (Gt + Sat), ^--J^[{y-1 + n)B,, + {a-3b + n)B^t + 2(a - 6 + n)Gl 
ri T 

^iii = ^{(« -b){y-l + n)B^t + 26(6 -2a- n)Bst - [2ab - (a - 6 + rO'lGj, 

b At 
Gt = xyBit — B2t, ai = a{ql), Ni = 



y — l + ri {z + xyriY 

z^, = z+- p, p = 2Ae J{t + z+)/V, Bu = Bi{ql, x[), i = l- 4. 



It is convenient to extract explicitly the contribution containing the infrared divergence. To do 
this we should add and subtract in the numerator of the integrand its value at z = z+. At this 
value we have: ri = 1, ai = a, Ni = Ar/y"^, x[ — x. So, the integral Iim can be written as 



1 dz r ^ , , oAt ^ , ,1 Vx ^ , p , doB , s 



Integral J^^n. Calculation of the integrand is performed in the approximation X2 = that 
corresponds to the emission of the collinear photon along the final-electron momentum. In this 
case the variables acquire the following values 

\-Z 2 _ ^2 



r2 = , ^2 = -r2QB: A 



l-z^'-' ■^--■^^ l-r2(l-2/) 

After integrating over the hard-photon angular variables, the integral /2m is represented as 
follows 

hm = 7^ r -^alN,J:,{z), (37) 

B JO Z-y — Z 
^2(2;) = T,2llRll + T,2ltRlt + ^2tt{Rtt — Rnn)i 



^2tt = hQlir^Gs+B^s), ^2it = -W^{(y-l+-)S4, + (a-36+-)S3,+2[l+(a-6)r2]G'J, 

V T r2 r2 

^211 = ^{(a -b)[l- r2(l - y)]B^s - 26[1 + (2a - 6)r2]S3, - [2abrl - (1 + ars - 6r2)']G'J, 

b At 

Gs = xyBu - 7^ zB2s, (^2 = a{ql), N2 = — ■ r^, -Bj, = Bi{ql, 4), i = 1 - 4. 

l-r2[l-y) [z + xyr2y 
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The contribution containing the infrared divergence is extracted exphcitly in a similar manner 
as it was done for the Iim integraL At the value ^ = 2;+ we have: r2 = 1, 0:2 = a, N2 = Ar/y"^, 
x'2 = X. So, the integral l2m is rewritten as 



1 /■^+ dz 



I- 2m 



At. 



Q% Jo z+ — y"^ J Try z+ dxdQ\ 



Vx , p , da 



The radiative corrections due to the virtual photon exchange and real soft- photon emission 
(with energy uj < Ae) can be related to the Born cross section as (note that the vacuum 
polarization effects are included in the Born cross- section through the dependence of the 
coupling constant a on the virtual-photon momentum) 

dxdQ\ dxdQ\ ' 

where the factor 5^^^ is IT 



«r/^ X, {^^f 3^ 1, 9,£l, TT^ 96', ^ ^ , Or 



here ei{e2) is the initial (final) electron energy, and 9 is the electron scattering angle in the 
coordinate frame /ci — ^2 + P = 0. The function / is defined as 

df 

f{x)= / -ln{l-t). 
Jo t 

The quantities ei{e2) and 6 can be expressed in terms of the invariant variables 

V{l-xy) V{l-z+) 2^ l-y-xyr 

£1 = — =, £2 = — =, cos - = — -. (41) 

2^V{t + z+) 2^V{t + z+) 2 {l-xy){l-z+) 

Then the radiative correction 6^^ is finally rewritten as 

^ = ^t-^"T~'^[ (i-xy)(i-.,) ]-^" (T3^) + (^-1)(3+2M (,_^^^^,_^^) ])}- 

(42) 

Integral Ir. When calculating this integral we use the results of Ref. [TH]. Besides the 
integrals calculated in that paper we need the following integrals 

F{z,r)xi, / 7— F(z,r)x?. (43) 



21TUJ J 2ttuj 

To calculate these integrals we write the hard-photon phase space in the form 

Qb dip , , 

= — ,Z —dzdr. (44) 

2TTUJ 2^y^ + 4a 2ti ^ ' 

Since, in our case, the function F does not depend on the ^9 variable, we can integrate over 
this variable. We do it in the above mentioned coordinate frame. The results are 



^1 = / 7^X1 



2-k'^ y^ + Aa 



i2-y)iy-c)-il-r)iy + 2a) , (45) 
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3if - 



1/2 + 4a 



(r — riY 



After simple calculations the integral Ir is (here we omit the contributions proportional to the 
Rin and i?^ components) 

1 r'"^^r+ rG'7™(^,r) G7*'^(^,ri)] 

1 — XV Jo Jr. 



r — ri 



1 — r 



1 — ri 



+ - 



1 



+ 



1 — 2;+ Jo Jr_ |r — r2 



+ 



+ 



(46) 



where 



\/7/2 + 4a Jo 

Li = /n 



1 — r 1 — r2 

Co7"(^, r) + iiCr(^, r) + i2Cr(^, r)] }, 



m?xy{T + 2;+) 
1 



, L2 = In 



m?xy{T + 2;+) 



(47) 



2xy{T + 2;+) 



2a;2/(r + z) + {z+ - z){y ± Jy^ + 4a) 



GT(;^,r) = ^(1 -r)5,Ar(^,r), G'r(^,r) = '^{l - r)B,B'r{z,r), 

with m,n = I, t, n. Note that the structure functions Bi are the functions of two independent 
variables = —rQ\ and ,x' = xyr / [z + ,T|/r). The expressions for the coefficients A™", i?™'™, 
C'fci", /c = 0, 1,2 are given in the Appendix A. The contributions proportional to the Ri^ and 
Rtn components are considered in a more detail in the Appendix B. 

Let us discuss briefly the singularities in the Ir integral. The value r — 1 corresponds to the 
real soft-photon emission (we have infrared divergence at this point), and the value r = ri(r2) 
corresponds to the emission of the collinear photon along the initial-(final-) electron momentum 
(the so-called collinear divergence). The singularity at the point 2; = 2;+ is the infrared one. 
The divergence at r = 1 point is unphysical one. It arises during the integration procedure due 
to the separation of the poles in the expression (xiX2)~^- It is necessary to extract explicitly 
the collinear and infrared divergences in the above formula. 

The integrand in the above expression can be written in the form which does not contain 
explicitly the infrared divergences if we add term (39) to it. To do this we use the following 
transformations 



In 



fi{x,y) iG{z,ri) r+ dr \G{z,r) G{z,ri) 



xy{T + z+) 



= P 



l-Ti 

dr 



+ 



L 



r — r,- 



1 -r 



1-r,; 



(48) 



r- [1 — r)\r — Tj 



G{z,r)-G{z,n) , i^l,2, 



where (pi{x,y) = (1 — xy)^, ip2{x,y) = (1 — z^)^, and the symbol P denotes the principal 
value of integral. The total radiative correction (which is the sum of the contribution due 
to the hard-photon emission and the contribution due to the real soft-photon emission and 
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virtual-photon contribution) to the part of the differential cross section caused by the tensor 
polarization of the target is written as 



da 



das 



dxdQg dxdQ^ 



+ 5'"\ 



(49) 



where 



-2/ 



5*°* = ^|3L + 2(L-l)/n 
27r 



da 



B 



ay 



l-xy){l-z^) 

dz 



4-^- 
3 



1 — xy) (1 — Zj^) -I dxdQ% xQ% Jo z^ — z 



8r 



afNiEi (z) + aiN2J:2{z) -a'y^E^ {z+) 



+ 



+ 



ay 



2xVQ^ 



^+ dz 



B 1=1 mn 

Q 



z+ — z 



Gr{z,n)-Gr{z+,1)-Gr{z,r2)+Gr{z+,1) + (50) 



y/y"^ + Aa Jo 



dz 



^+ , a 



dr—Bi 



Cr{z,r)+iiCr{z,r)+i^Cr{z,T)\+Rr}. 



The term EJ^^ has different form depending on the integration region of the variable r. For the 
regions r_ < r < ri and r2 < r < r_|_ the function i?™" has the following form (in these regions 
r 7^ 1 and therefore the divergence at the point r = 1 is absent) 



1 / dz -[Gr{z,r)-Gr{z,n) 

1 — xyJo Jr- (1 — r)|r — ri| L 



+ 



(51) 



+ 



dz 



r+ 



dr 



1 — z^ Jo 
For the region ri < r < r2 we have 



[1 — r) |r — 



Gr{z,r)-Gr{z,r2) 



Rr - dzln(^—^) {gTiz, 1) - f^iz, 1) + ^ 
Jo — 1^ z+ — z 



g^-{z,l)-gT{z,n) + 



1 ^ /'''+ dv r 

+ friz, 1) - fTiZ: r2)] } + l dz I r^igTiZ: t) - gT{z, 1) - 



(52) 



where we introduce the following notation 

Griz, r) = gZiz, r) + A,g^-{z, r), G^iz, r) = f^iz, r) + A2/r 0, 



F"^"(z,r) 



1 



r — ri 



9TM-gT{z.n) , F^-{z,r) 



1 



r — r2 



fT{z.r)-fT{z.r2) , (53) 



Ai = (1 — a;|/)r — a — 6 — 2;, A2 = (l — y + a;y)r + ^ — 1. 
In obtaining the above formula we use the relation 



i - r 



1 — 



>f(r) - <^^(i)] + ^{l)ln( -) 



(54) 



At last, let us consider the part of the integral / which is caused by the Rin and Rtn 
components of the deuteron quadrupole-polarization tensor. As stated above, these components 
do not contribute to the cross section treated in the Born approximation. If these terms are 
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integrated over the whole region of the tp variable, then these integrals are equal to zero as 
well (this result is due to the fact that only one plane is remained after such integration). We 
discuss this problem in more detail in Appendix B. 

Note that the integration limits in formula (50) over the variable z are given somewhat 
schematically. This integral contains two contributions (we neglect here the contribution of 
the radiative tail from the quasi-elastic scattering). One of them is the so-called inelastic 
contribution and the integration region for it over the variables r and z is presented in Fig. 1 
by the dashed triangle. The integration over z variable for this contribution must be carried 
out from Zmin = {M'^h ~ M)/V to 2;+, where Mth is the inelastic threshold {Mth = M + m^r). 
The second contribution, related to the radiative tail of the elastic peak, is given by the interval 
z = 0, r_(0) < r < r+(0). 




Figure 1: The integration domain with respect to the r and z variables. 

The contribution of the elastic radiative tail (i.e. the inclusion of radiative corrections to 
the elastic ed-scattering) to the total radiative correction 5*°* can be obtained from formula 
(50) by simple substitution in the hadronic tensor 

BM\x')^--^5{l-x')Bf\ i = l-4, (55) 

where are expressed in terms of the deuteron electromagnetic form factors as 

B^^ = m'Gl,, = -2r^\^Gl, + i^(Gc + ^Gq + t^Gm)], (56) 

L -\- Tj o 
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Here Gc, Gm and Gq are the deuteron charge monopole, magnetic dipole and quadrupole 
form factors, respectively. These form factors have the following normalizations: 

Gc(0) = 1, Ga/(0) = (M/m„)/i,, G'q(O) = M^Q,, 

where m„ is the nucleon mass, HdiQd) is the deuteron magnetic (quadrupole) moment and their 
values are: /x^ = 0.857, Qd = 0.2859/m^. After substitution of i^j-*^'^ into formula (50) we have 
to do a trivial integration over z variable using a (5-function 5(1 — x') = xyr6{z). 



4 Numerical estimations 



We calculate the radiative corrections for the kinematical conditions of the HERMES ex- 
periment |T0|. The energy of the positron beam is 27.6 GeV. The HERMES has provided the 
first direct measurement of the structure function bi in the kinematic range 0.002 < x < 0.85 
and O.lGeV^ < < 20GeV^. A cylindrical target cell confines the polarized gas along the 
positron beam line, where a longitudinal magnetic field provides the quantization axis for the 
nuclear spin. The corresponding tensor atomic polarization is T=0.83 (for the definition of this 
quantity see the Appendix C). 

The analysis of the experimental data was performed in the approximation 63 = 64 = 0. In 
further numerical estimation we also neglect these functions. 

The deuteron spin-dependent structure function hi is extracted from the measured tensor 
asymmetry A^z via the relation ^U] 

b - \ (^7) 
- ~2^^^ 2x{l + R)' 

where the deuteron spin-independent structure function has been expressed in terms of the 
ratio R = ai/ar = F^^l + AM'^x'^ /Q'^)/2xFf-l [20], 7 is a kinematic factor (7^ = AM"^ x'^ / Q'^) . 
Here (Jt{o'l) is the cross section for the absorption of the transversely (longitudinally) polar- 
ized virtual photons by the unpolarized target. The Born cross section of the deep-inelastic 
scattering of unpolarized electron beam by unpolarized target has the form 

- -[{l-y-xyT)F^{x,Q^) + xy'Ft{x,Q% (58) 



The structure functions Ff^ are related to the structure functions Wi^2 (introduced in the 
formula (7)) by the following way: Wi = 2Ff, W2 = 4(r/j/)F^. The deuteron spin-independent 
structure function F^ = F|(l + F^/F^)/2 is calculated using parameterizations for the proton 
structure functions s21^ and the ratio F^ IF2 |221- The deuteron spin-dependent structure 
function 62 has also been extracted from the experiment using the Callan-Gross relation 



62 = 2x^^6i. (59) 



1 + i? 

1 + 7 

According to the preliminary results of the HERMES experiment the tensor asymmetry can 
be parametrized as 



A,, = -1.56 ■ 10^2(1 - 1.74x - 1.45v^). (60) 
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Figure 2: The spin-dependent part of the cross section calculated for the kinematical conditions 
of the HERMES experiment ^U]. The solid line is the Born approximation, and the dotted 
line corresponds to the inclusion of the radiative corrections. The values are the following: 
a - O.lGel^^ h - IGeV^, c - AGeV^, d - lOGeV^. 

The influence of the radiative correction on the spin-dependent part of the Born cross- 
section is shown in Fig. 2 as a function of the variable x for various values. The inclusion 
of the radiative correction shifts the zero value of bi and 62 to the smaller x-value region (see 
Fig. 2 c and d). In the range of low x (x ~ 10^'^ — 10^^) the value of radiative correction 
changes from 10% to 30% as compared with the Born contribution. It is the region that of the 
utmost importance for 61-measurements. According to the theoretical predictions 
the structure function bi increases very rapidly in this region and this fact was confirmed in 
the HERMES experiment [Til] . 

From our estimation we conclude that the radiative corrections to process (1) are not small, 
especially for low x region, and they have to be taken into account in the data analysis. 
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5 Appendix A 



In this Appendix wc present the formulae for the coefficients A^"', B^"', and C^'^, {m, n = 
l,t, i = 1 — 4, j = 0,1,2) that determining the cross-section of the hard-photon emission 
process (see formula (50)). 

The coefficients determining the contribution proportional to the Bi structure function are: 

= -^[(f - Ai)2 - 2a(b + Ai)], = ^{[(2a -b)r + l + A2)]' - ar(2 + 3ar)}, 

C = -^{(^-Ai)' + a[3a(l + r^)-2(6 + Ai)]}, Ci\ = -67V(c + 2a), = -6A:^, 

4* = M^QlC^b + AO(f - AO, = J^QU^2 - 2br)[{a -b)r + l + A,], (A.l) 
C^*^ = -4:n2Ql[a{l+r^){y+2a)-2bf-Ai{c+2a-2b)], C[\ = -4^2 [(|/+4a)(c+2a)-2a(f+26)' 

C^*i = -8n2^(y + 2a), A'^ ^ -^^[b' + (b + A^f], Bf = ^^{2b'r' - 2brA, + Al), 



ritt _ 



2b 



2 N 

(l+r')(y2+4a-2a&) + (26+Ai)2 + A?], C« = — [6(l+y+2a-r) + (l+a)Ai 



+ 2a{2 - b) 



The coefficients determining the contribution proportional to the B2 structure function are: 



.11 ^3 



6(1 + r^) + (1 - r + ry) Ai (f - Ai)^ - 2a{b + Ai 



r 
A\/ 



6(1 + r^) - A2(f - 2a)] {[(2a - 6)r + 1 + Aa]^ - ar(2 + 3ar)}, 



= -^{(7-3y)c^ + 3a(5-|/ + r)c + 3a2(3 + r2) -ar[5 + 3(a + 6)^]}, 
Cf2 = -3t- [4(a + c) - cy] , C^^ = "6^^' 
4* = -2^(2& + Ai)(f - Ai) [6(1 + r^) + (1 - r + ry)Ai 
B'^ = -2^^(A2 - 26r)[l + A2 + (a - b)r] [6(1 + r^) + (a + 6 - r)A2 



(A.2) 
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C'o\ = ^{2ac[yf+{3b+a){l+r)-y-8a\- 



c 



2a+{2-y){y+4a) 



+2a 



2a{b-a+r) + {y+2a){r- 



-a{l+r'^)+r{a+bY)]} , Cj^ = -4n2- [2a(l-3r+2f+4c)+c?/(2+6-a) 



, eg = -8n2^{y+2a), 



^« = ^nsZ,[2b^ + Ai(26 + Ai)], S« = -^^n3Z2[26V + A2(A2 - 2br)], 
b b 



xy 



VN 



C** = -^[-2c^[a+{l + a){2-y)\ + c[(3 - 2?/ + + 6^)(f - 2a) +4(a6 + 6- a^) +4r(a - 6^) 



-2a 



(r - a)^ + 6^] + (1 + 2a - 26 + + 6^) [r - a(l + r^) + (a + }, 



2cd^ 



{c[l + 7a(l + a) - 6(1 + 6) + (a + 6) (a^ + 6^ 

iVr 



+ 4a 



:a-6)(l-r) + a2 + 62-r]}, 



L/99 



1/^ + 2a(2 - b) 



The coefficients determining tlie contribution proportional to tlie B3 structure function are: 

A^^ = n3-{(a + f)[2Zi + rAi(2a + r-Ai)] - Ai p(r -Ai) + 2(6 + Ai) +r(a + 6) (a + r -Ai)] }, 

^" = -713^(2^2 [1 + (2a - 6)r + A2] + 3aA2[(6 - a)r - 1 - A2]}, 



'3 - -^3 

VN 



Cos = ^ [c(6a - 16 + 9y) + 6a(y - 3 - r)] , Cfa = -3t7V(2 - y), C| = 0, 



2 

cQ' 



(A.3) 



A^* = J^{2Zi(36 - a - r) + Ai [4r(l + 6^ + 3a6) - 2a(l + r^) + c(ar - 3 + 56r)] }, 

IVJ. (Jj 



a)r — 1 



+ A2 (1 + ar)(a - 66) - (a + 36)(r2 + A2) + r{b^ - 1) + 6+ 



+A2(3r - 26)]}, C^3 = n2y{4a[26f - + 4(6^ - a)] - c[3|/(2 - y) + 8a(l + a + 26)]}, 

Cf3 = -2iV^(2 - y)(y + 2a), C| = 0, 
4* = n3a;7/^{(Ai - 26)[6(1 + r^) + {l-r + ry)A,] + 6Ai[l + r(6 - a + Ai)]}, 
Bf = nsxy^{i2br - A2MI + r^) + (1 - r - y)A2] + 6A2[r(6 - a + r) - A2)]}, 

C« = ^{36-a-(a2+62)(2+a+6)+r[yV2y(26-a)+2a(3-a)]+c[y(l+y+3a)-4(l+a)-2a6]}, 



^13 



VN 
2d? 



{2-y)[y' + 2a{2-b)], C«3 = 0. 



Tlie coefficients determining the contribution proportional to the structure function are: 

2 2 
A^l = 713-1(6 - a)(l + r') + Ai [l + r(2a - 6)] |, B^^ = 713- f(a - 6)(1 + r') + A2(a + f) 



r AT 



4* = ^3^{26(1 + r') + Ai [1 - r(36 - a)] }, = 713 
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l + 3(6-a)], Cf4 = C^4 = 0, 

-26(l + r^) + A2(f-26) 



Md 



At = -xyrc^Aina, Bf = -xyc^AiUs, C^l = 



1 + 4a6 - (a - 6)^] , C[\ = = 0, 
cNV 



(A.4) 



(y + 2a), C« =C« =0, 



here we use the following notation 



C — Z + xyr, r — a — b + r, rii 



~2 1-r 



VQl, 712 



NV 
2Md 



= hQ%^ Ai = (1 - xy)r -a-b- z, As = (!-?/ + xy)r + z - 1, N = 
Zi = 6(1 + r^) + Ai(l - r + yr), = b{\ + r^) + A2(l -y-r). 



N 1/2 
2^1 -r' 

4r 



6 Appendix B 

In this Appendix we consider the part of the integral / which is caused by the and 
Rtn components of the deuteron quadrupole-polarization tensor (these components do not 
contribute to the differential cross section treated in the Born approximation). Let us define 
the integral which is caused by the component 



27ra; 



T,in{z,r,ip)Rin, 



(B.l) 



with 



Piin = --^{cgxy{l + r'^)Bi + g\c{ 
1 — r L 



Xl X2 

1 - r(l - y)) + a(l + r') - 4/r]52 + [2a - fr+ 



+ -(3(1 - r + yr) + 2ar) 



S3 + ^c[l + (a-6)r]S4}, (B.2) 



Piin = -^^{-xy(l + r2)(c + 2ar)Si + -[-ar(a(l + r^) - 4/r) + ^(1 - y - r)+ 
1 — r c z 

+|(4/r - a(l + + 2yr - 2r))]52 + [r(/ - 2ar) - |(2a + 3(r + y - \))\B^ - |(a - & + r)54}, 
C/o,n = 2g(cBr + + 2t(2 - y)(B2 + S3), C/i^n = y (^i + ^^2), 



and c = z + xyr, / = 1 + (1 — y)^, (? = 1 + 2a/c, n ■ q = Sjj^^q^. 
The second integral, caused by the Rtn component, is defined as 

f d^k 

where the integrand is 

a^(g^) 2VN (Pun P2tn^jj ^ 



(B.3) 
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Pitn = Ql{fxy{l + r^)B^ - ^ [a{l + r') - 4r(/ + 4y)]B2+ 
+/[1 + r{y - 1)]{B2 + Bs) + hrc{B^ + B^) + 26r(l - y)B^}, 
Pnn = -Ql{xyg{l + r')B, + ^[a(l + r') - 4fr]B2+ 
+ g{r-l+y) {B2 + S3) - 60(^3 + S4) - 2fer(l - y)Bs}, (B.4) 
Uotn = (r - 1) \-2xyf{B, + -^2) + (2 - |/)(2a + y){B2 + ^a)], 



C/un = :^{(r-l)(2a + y)Gi-/G2-|^[M(2a + y)-2rd](S2 + S3)-^[c-2a(r-2)]G2}, 

Gi = 3Si + —B2 + -^{B2 + 2B3 + S4), G2 = -Si - -^{B2 + 2S3 + S4), 

and (i^ = bQ'^^, f = b — a — z + r{l — xy), g = z — l + r{a — b + xy). 

As before, we calculate the above integrals in c.m.s. of the hard photon and undetected 

— * — * — * — * 

hadron system: ki — k2+p — 0. The electron momenta ki and k2 define the xz plane, z axis is 
directed along the deuteron momentum p. Then the hard-photon momentum k is determined 
by the azimuthal {(f) and polar {9) angles, and the phase space of the hard photon can be 
written as 

pzdr, (B.5) 
27ra; 2y/y^ + Aa 27r ' ^ ^ 

where cu is the hard-photon energy. 

The quantity n • q can be written in this coordinate system as n • q — nsinif, where ri is a 
factor independent on (p. Then the integration over the f variable in region (0, 27r) leads to the 
following result: = Ifn = 0. So, the and Rfn components of the deuteron quadrupole- 
polarization tensor do not contribute to the differential cross section of deep-inelastic scattering 
of unpolarized electron beam by the tensor polarized target. This result is due to the fact that 
only the scattered-electron variables are measured (it corresponds to the HERA experimental 
conditions, for example). 

If the hard-photon is detected then the and survive and the expressions for E;„ and 
'Etn have to be taken into account 



7 Appendix C 



In this Appendix we give some formulae describing the polarization state of the deuteron 
target for different cases. For the case of arbitrary polarization of the target it is described by 
the general spin-density matrix (in general case it is defined by 8 parameters) which in the 
coordinate representation has the form 

Pfj-v — ^(fi'/ii/ -^2 ) 2M^^'^''^'^^^'' QiJ-fi QiJ-v ~ Qvixi Qnn ~ ) PixQfiv — , (c.i) 

where p^ is the deuteron 4- momentum, and Q^i, are the deuteron polarization 4- vector and 
the deuteron quadrupole-polarization tensor. 
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In the deuteron rest frame the above formula is written as 

1 i 

Pij = 3% + ^^ijkSk + Qij, ij = X, y, z. (C.2) 

This spin-density matrix can be written in the hehcity representation using the following rela- 
tion 

Pxx' = Pije^^*ef'\ A, A' = +, -, 0, (C.3) 

where e^^^ are the deuteron spin functions which have the deuteron spin projection A on to the 
quantization axis (z axis). They are 

eW = T-^(l,±z,0), e(°) = (0,0,1). (C.4) 

The elements of the spin-density matrix in the helicity representation are related to the ones 
in the coordinate representation by such a way 

111 1 1 

P±± ^ 3 T 2"^^ ~ 2^zz, Poo = 3 + Qzz, P+- — —-^{Qxx — Qyy) + ^Qxy, (C-5) 
11 11 

P+0 = -:^^{sx-iSy)---j={Q^:,-iQy^), p^Q = -^^^{sx+iSy)+--^{Q^^+iQy^), px\> = (pa'a)*- 

To obtain these relations we use Qxx + Qyy + Qzz = 0. 

The polarized deuteron target which is described by the population numbers n+, n_ and Uq 
is often used in the spin experiments. Here n+, n_ and Uq are the fractions of the atoms with the 
nuclear spin projection on to the quantization axis m = +1, m = —1 and m = 0, respectively. 
If the spin-density matrix is normalizd to 1, i.e. Spp = 1, then we have n+ + n_ + no = 1. 
Thus, the polarization state of the deuteron target is defined in this case by two parameters: 
the so-called V (vector) and T (tensor) polarizations 

y = n+ - n_, T = 1 - 3no. (C.6) 

Using the definitions for the quantities n±,o 

(±)* (±) (0)* (0) 

n± = pijel ' e) \ no = pije\ e) \ (C.7) 

we have the following relation between V and T parameters and parameters of the spin-density 
matrix in the coordinate representation (in the case when the quantization axis is directed 
along the z axis) 

1 111 
^0 = 3 + Qzz, n± = - T -Sz - -Qzz, (C.8) 

or 

T = -3Q,z, V = -s,. (C.9) 
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